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ABSTRACT 

Network flows is an important subfield of Operations Research and is 
well-known for its diverse applications and highly efficient algorithms. 
Applications of network flow models have been identified in a variety of 
situations including applied mathematics, communication systems, defense, 
manufacturing and production, scheduling, and public policy. In the recent 
past, researchers have made efforts to identify, compile and summarize 
those applications which can be transformed into fundamental network flow 
problems such as shortest paths, maximum flows, minimum cost flows, 
assignments and minimum spanning trees, etc. This thesis addresses a 
similar issue and attempts to identify the applications of network flow 
models in Lagrangian Relaxation. 

Lagrangian relaxation is an important technique in combinatorial 
optimization and is extensively used in designing branch and bound 
algorithms for hard (i.e., NP-complete) problems. The relaxed problems 
obtained by applying the Lagrangian relaxation technique are often network 
flow problems and, therefore, can be very efficiently solved by network 
flow algorithms. 

In this work, we make an attempt to identify and summarize such 
applications whose lagrangian relaxation results in network flow 
subproblems. We report over 40 applications which have been selected 
through an extensive research of the available literature. These also 
include applications which we have formulated ourselves. We hope that the 
understanding of the models described here will allow researchers to 
identify more and more hard optimization problems whose lower bounds can be 
obtained more efficiently by using network flow algorithms in Lagrangian 


relaxation. 



IV 




CONTENTS 

Page 

Chapter 

i 

INTRODUCTION 



i.i 

Introduction 

1 


1.2 

Notation and Definition 

4 


1.3 

Network flow problems 

6 


1.4 

Overview of the thesis 

11 

Chapter 

2 

LAGRANGIAN RELAXATION 



2. 1 

Introduction 

13 


2.2 

Literature Survey 

13 


2.3 

Lagrangian Relaxation Technique 

14 

Chapter 

3 

APPLICATIONS OF SHORTEST PATH PROBLEMS 

21 

Chapter 

4 

APPLICATIONS OF MINIMUM SPANNING TREE PROBLEMS 

37 

Chapter 

5 

ADDITIONAL APPLICATIONS 

49 



REFERENCES 

70 



1 


CHAPTER 1 

INTRODUCTION 


1 . 1 INTRODUCTION 

Network flows is an important area in linear programming and is known 
for its novel theoretical properties. For most of the network flow 
problems, very efficient algorithms have been developed by exploiting the 
structure of the constraint matrix. However, many real life and practical 
problems can be viewed as easy network flow problems complicated by a 
relatively small set of side constraints. These problems can not be 
modeled as pure network flow problems and are difficult to solve as direct 
algorithms are not available for such constrained network flow problems. 

Lagrangian relaxation is a flexible solution strategy that permits to 
exploit underlying structure in any optimization problem by relaxing 
complicated constraints. This approach permits to "pull apart" constrained 
network flow problems by removing complicating non-network constraints and 
instead place them in the objective function with associated Lagrangian 
multipliers. Dualizing these complicating side constraints using the 
Lagrangian relaxation technique produces Lagrangian subproblem (or 
subproblems) that will be a pure network problem in nature. This 
Lagrangian subproblem can be solved as a pure network flow problem and its 
optimal value will be a lower bound (for minimization problem) on the 
optimal value of the original constrained network flow problem. Many 
practitioners identified various applications of such hard combinatorial 
optimization problems which are having two sets of constraints, one set of 
“network" constraints and the other of non-network side constraints. These 
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applications are few in number and appear to be scattered in the 
literature. 

Though Fisher [1981] and Ahuja, Orlin and Magnanti [1993] gave a 
compilation of few such applications, yet there is no single paper, or any 
other reference that summarizes all such real life applications at one 
place. This report fills precisely this need by describing, summarizing or 
simply referencing over 40 such combinatorial optimization problems in 
which by relaxing complicating non-network constraints, we can get pure 
network flow Lagrangian subproblems. Furthermore, we also tried to 
identify and add some applications of similar nature from our side in this 
report. The network flow problems, which arise as Lagrangian subproblems 
in such relaxed hard combinatorial optimization problems, are: 

(i) The shortest path problem. 

(ii) The minimum spanning tree problem. 

(iii) The minimum cost flow problem. 

(iv) The assignment problem. 

(v) Convex cost flow problem. 

(vi) Generalized flow problem. 

(vii) Multicommodity flow problem. 

Identifying such constrained network flow application or combinatorial 
optimization problems with two sets of constraints, proved to be a 
difficult task. Since there existed no reference on such applications, we 
almost had to survey the entire combinatorial optimization/network flows 
area. 

In the first phase, we went over the papers cited in the Integer 
Programming Bibliographies (Lasting [1976], Hausman [1978] and Von Randow 
[1982, 1985]) under the subjects lying within relaxation, network flows and 



3 


integer programming. We also went through the international abstracts of 
operation research in above areas. We then scanned these papers to 
identify such applications. We also scanned the papers given in the 
reference list of above identified papers. In the second phase, we went 
through many journals in management science, operations research, 
mathematics and communications where such articles have appeared in the 
past. We also thought some additional applications from our side in this 
phase. 

During our research, the applications we found, could be categorized 
into two broad classes: 

(i) Applications whose Lagrangian relaxation resulted into one single 
pure network flow problem. 

(ii) Applications whose relaxation resulted into more than one, different 
type of pure network flow subproblems. 

We have included in our report most of the such applications that were 
interesting, simple to describe and throw insight on modeling and relaxing 
techniques. We have kept the readability and the importance of an 
application in mind while deciding its inclusion in the text. There were 
another category of applications that were nice but either too related to 
the applications that we described completely or rather simple to be 
discussed fully. For such applications, we only describe the problem. To 
summarize, the applications we identified are covered in this report in the 
following three manners: 

(i) The problem is described and its formulation is shown. By using 
Lagrangian relaxation, resulting pure network flow subproblem (or 
subproblems) is shown for such application. 

(ii) The problem is described and its equivalence with previously 
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described problem has been proved. 

(iii) The problem is simply stated with no description of the model given, 
because of the complexity of the problem. 

This report serves the following purposes. For the first time it 
collects together all such real life applications which are network flow 
problems with additional side constraints. These application will throw 
insight on the use of Lagrangian relaxation technique to solve such hard 
problems. These applications are of great use to teachers who may include 
some of these in their courses in combinatorial optimization, Operations 
Research and Network Flows. We believe that the understanding of the 
models described here will allow the researchers to bring more and more 
hard optimization problems those can be solved as network flow problems 
using Lagrangian relaxation. 

1.2 NOTATION AND DEFINITION 

We now collect together several basic definitions and describe some 

notation. We consider a directed graph G = (N,A) consisting of a 'set, N, 

of nodes, and a set, A, of arcs whose elements are ordered pairs of 

distinct nodes. A directed network is a directed graph with numerical 

values attached to its nodes and/or arcs. We let n = j N | and m = | A | . We 

associate with each arc (i, j) € A, a cost c. . and a capacity u. . £ 0. 

J ^ J 

Frequently, we distinguish two special nodes in a graph: the source s and 
sink t. 

An arc (i,j) has two end points, i and j. We refer to node i as the 
tail and node j as the head of arc ( I , j ) , and the arc (i,j) is incident to 
nodes i and j. The arc (i,j) is an outgoing arc of node i and an incoming 
arc of node j. The arc adjacency list of node i, A(i), is defined as the 
set of arcs emanating from node i, i.e., A(i) = I(i»j) € A:j € N>. The 
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node-arc incidence matrix representation, or simply the incidence matrix 
representation, represents a network as the constraint matrix of the 
problem. This representation stores the network as an nxm matrix [N which 
contains one row for each node of the network and one column for each arc. 
The column corresponding to arc (i,j) has only two nonzero elements: it has 
a +1 in the row corresponding to row i and a -1 in the row corresponding to 
row j. The degree of a node is the number of incoming and outgoing arcs 
incident to that node. 

A directed path in G = (N,A) is a sequence of distinct nodes and arcs 
ij, (ij, i 2 ), i 2 , (i 2> i 3 ), i 3> . . - , (i r _ 1 , i r ), i r satisfying the property that 
(i k ,i k+1 ) e A for each k = l,...,r-l. An undirected path is defined 
similarly except that for any two consecutive nodes i^ and i^ +1 on the 
path, the path contains either the arc (i , i^ ) or the arc ^k+l ’ 1 k^ ‘ A 
directed cycle is a directed path together with the arc (i .i^) and an 
undirected cycle is an undirected path together with the arc (i^.i^) or 

(i r i r ) - 

We shall often use the terminology path to designate either a 
directed or an undirected path, whichever is appropriate from context. For 
simplicity of notation, we shall often refer to a path as a sequence of 
nodes i - i - ... - i, when its arcs are apparent from the problem 

JL 

context. Alternatively, we shall sometimes refer to a path as a set of 
(sequence of) arcs without any mention of the nodes. We shall use similar 
conventions for representing cycles. 

A graph G = (N,A) is called a bipartite graph if its node set N can be 
partitioned into two subsets and N 2 so that for each arc (i,j) in A, i e 
N 1 and j 6^. 

A graph G' = (N',A' ) is a subgraph of G = (N,A) if N' £ N and A' £ A. 
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A graph G 7 = (N 7 ,A 7 ) is a spanning subgraph of G = (N,A) if N 7 = N and A 7 Q 
A. 

Two nodes i and j are said to be connected if the graph contains at 
least one undirected path from i to j. A graph is said to be connected if 
all pairs of its nodes are connected; otherwise, it is disconnected. The 
connected subgraphs of a graph are called its components. We always assume 
that the graph G is connected and hence m £ n-1. 

1.3 NETWORK FLOW PROBLEMS 

We shall now briefly describe the network flow problems which arise as 
Lagrangian subproblems in the relaxed problems considered in this thesis. 

Shortest Path Problem 

The shortest path problem is to determine directed paths of smallest 
cost from a given node s to all other nodes in the network. The problem 
can be stated as the following linear programming problem: 


minimize } C. . X. . 

L xj ij 

( i , j ) eA 


subject to 


Y X. . - , 

L ij L Jl 

{ j: (i, j)eA> { j: ( j,i)eA} 


l *- ■ < 


f(n-l) for i = s, 

-1 for all i € N -{s>, 


X. . i 0, for all (i,j) € A. 
ij 


Some of the simplest applications of the shortest path problem are to 
determine a path between two specified nodes of a network that has a 
minimum length, or a path that takes least time to traverse, or a path that 
has the maximum reliability. 


Assignment Problem 
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The data of the assignment problem consists of a set, N^, say of 
persons and a set N 2> say of objects, satisfying = |N 2 |, a collection 
of node pairs A S x representing possible person-to-object 
assignments and a cost C_ associated with each element (i,j) in A. The 
objective is to assign each person to exactly one object in a way that 
minimizes the total cost of assignment. This problem can be stated as the 
following linear program: 


I C,,X,, 


minimize ^ ij ij 

(i, j) eA 


subject to 


V X. . = 1 for all i e N. , 

L ij 1 

{j: (i, j) €A > 

£ x ij = i for all j e N 2 , 
{j: (i, J^ eA} 


X.j s 0, for all (i,j) e A. 


Sometimes it may be unba l anced problem that can be transformed to a 
balanced assignment problem by introducing some dummy nodes. Sometimes, we 
have a utility u.j assigned with each arc (i,j) €A and our objective is to 
obtain an assignment that maximizes the total utility £ u. , X_. We 


( i , j ) eA 


ij iJ 


can 


transform this problem to the previous version by defining Ch = - u.^ 


and minimizing the objective X u X.. • 

( i , j ) eA 


The applications of the assignment problem include assigning people to 
projects jobs to machines, tenants to apartments, swimmer to events in a 
swimming meet, and medical school graduates to available internship. 
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Minimum Spanning Tree Problem 

A spanning tree is a connected acyclic graph that touches all the 
nodes of an undirected network. The cost of a spanning tree is the sum of 
the costs of its arcs. In minimum spanning tree problem, we wish to 
identify a spanning tree of minimum cost. 

Minimum Cost Flow Problem 

We wish to determine a least cost shipment of a commodity through a 
network that will satisfy demands at certain nodes from available supplies 
at other nodes. The problem can be stated as the following linear 
programming problem: 

minimize ) C. . X. . 

L ij ij 

( i , j ) eA 

subject to 


Y x. . 

L xj 


L ij u j 

{j: (i, j)eA} {j: (j, i)eA> 


Y X.. = b(i) , 

L ji 


1. . * X. . s u. . . 

ij iJ ij 

In this problem, C_ represents the cost per unit flow on arc (i.j), 

1. . represents the lower bound on the arc flow and u. represents upper 

ij 1 * * * * * * * * X J 

bound on arc flow. The number b(i) represents the supply/demand of a node 

i e N. If b(i) > 0, then node i is a supply node; if b(i) < 0 then node i 


is a demand node; and if b(i) = 0 then it is a transshipment node. 


The applications of minimum cost flow problem include : . the 
distribution of a product from manufacturing plants to warehouses, or from 
warehouses to retailers; the flow of raw material and intermediate goods 

through the various machining stations in a production line; the routing of 
automobiles through an urban street network; and the routing of calls 
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through the telephone system. 

Convex Cost Flow Problem 

In the minimum cost flow problem, we assume that the cost of the flow 
on any arc varies linearly with the amount of flow. Convex cost flow 
problems have a more general cost structure: the cost is a convex 
function of the amount of flow. The problem can be formulated as follows : 

minimize T C (X ) 

( i , j ) g A 1J 1J 

subject to 

£ X. . - £ X.. = b(i ) , for all i € N. 

(j:U,j)€A) 1J ( j: ( j , i )eA) J1 

0 s X. . < u. for all (i,j) e A, 
ij lj 

X is integer for all (i,j) e A. 

We define this model on a directed network G = (N,A) with a capacity 

u and a convex cost function C. .(X. .) associated with every arc (i, j) 6 A. 
X J ij ij 

As always, we associate a number b(i) with each node i € N specifying the 
node’s supply or demand, depending upon whether b(i) > 0 or b(i) < 0. 

Flow costs vary in a convex manner in numerous problem settings 
including (i) power loss in a electrical network due to resistance; (ii) 
congestion costs in a city transportation network; and (iii) expansion 
costs of a communication network. 

Generalized Flow Problem 

In the minimum cost flow problem, arc conserve flows, i.e., the flow 

entering an arc equals the flow leaving the arc. In generalized flow 

problems, arc might "consume" or "generate" flow. If X„ units of flow 

enter an arc (i,j), then u. .X. . units arrive at node j; u. . is a positive 

ij IJ iJ 

multiplier associated with the arc. If 0 < < 1, then the arc is lossy, 
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and if 1 < P iJ < «, then the arc is gainy. The problem can be formulated 
as follows : 

minimize V C. . X. . 

(i.J)€A 1J 1J 

subject to 

I X. . - I p.. X., = b(i) for i e N, 

<j:(i,j)€A} 1J ij: ( j, i )eA> J1 J 

0 s X. . s u, ,, (i» j) e A. 
ij 1J 

Generalized network flow problems arise in several appli ca ti° n 
contexts: for example (i) power transmission through electric lines when 
power is lost with distance; (ii) flow of water through pipelines or 
canals that lose water due to seepage or evaportion; (iii) cash management 
scenarios in which arcs represent investment opportunities and multipliers 
represent appreciation or depreciation of an investment’s value. 

Multicommodity Flow Problem 

The multicommodity problem arises when several commodities use the 
same underlying network. The commodities may either be differentiated by 
their physical characteristics or by simply their origin-destination pairs. 
Different commodities have different origins and destinations, and 
commodities have separate mass balance constraints at each node. However, 
the sharing of common arc capacities binds the different commodities 
together. In fact, in multicommodity problems we deal with the allocation 
of the capacity of each arc to the individual commodities in a way 
minimize the overall flow cost. 

Let X^. denote the flow of commodity k on arc (i,j)> and let 

X J 

( ’ * Th.Gi 

denote per unit cost of transportation of commodity k on arc U,J1- 
problem can be formulated as follows : 
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K 


k k 
C K . XV. 


minimize V V 

k=l ( i , j )eA 1J 1J 


subject to 


y — V x^ ^ 

{j:(i,j)eA> ^ {j:(j,i)eA> ^ * 





for all 

i e N, k = 1 K, 


K 

I 

k=l 

X . . - 
ij 

u. . 
ij 

for all 

( i , j ) e A, 


0 * 

x k . s 

k 

u. . 
ij 

for all 

( i , j ) e A and k = 1,2,. 

. . ,K. 

The total flow of 

all 

the 

commodities 

on each arc (i,j) has 

been 


restricted to maximum arc capacity u„. Furthermore, individual flow 
k 

bounds on the flow of commodity k on arc (i,j) has been imposed. 

Multicommodity flow problems arise in many practical situations 
including (i) the transportation of passengers from different origins to 
different destinations within a city; (ii) the worldwide shipments of 
different varieties of grains from countries that produce grain to those 
which consume it. 


1.4 OVERVIEW OF THE THESIS 

In Chapter 2, we describe the Lagrangian Relaxation technique and 
related properties and theorems. In Chapter 3, we describe applications 
in which relaxed Lagrangian subproblems are shortest path problems. 
Chapter 4 includes applications in which relaxed Lagrangian subproblems are 
either undirected minimum spanning tree problems or directed minimum 
spanning tree problems. In Chapter 5, we describe applications whose 
Lagrangian relaxation result in more than one Lagrangian network flow 
subproblems. In addition to these, we also report the additional 
applications in which Lagrangian subproblems will be network flow problems 
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other than shortest path and minimum spanning tree problems. 

We have mentioned references against each application. We have not 
mentioned any reference against applications which we have formulated 
ourselves. 
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CHAPTER 2 

LAGRANGIAN RELAXATION 


2.1 INTRODUCTION 

Combinatorial optimization problems come in two varieties. There is a 
small number of easy problems which can be solved in time bounded by a 
polynomial in the input length and an all-too-large class of "Hard 
Problems" for which all known algorithm require exponential time in the 
worst case. Many hard problems can be viewed as easy problems complicated 
by a relatively small set of constraints. 

Lagrangian relaxation is a solution strategy that dualizes side 
constraints and produce a Lagrangian subproblem that is easy to solve and 
whose optimal value is a lower bound (for minimization problem) on the 
optimal value of the objective function of the original problem. The 
Lagrangian relaxation can thus be used in place of a linear programming 
relaxation to provide bounds in a branch and bound algorithm. 

2.2 LITERATURE SURVEY 

The Lagrangian multiplier technique of nonlinear optimization dates to 
the 1700s and was originated by the famous mathematician Lagrange, from 
whom the technique is named. There were a number of forays prior to 1970 
into the use of Lagrangian methods in discrete optimization, including the 
Everett's proposal for "generalizing" Lagrange multipliers (1963). 
However, the "birth" of the Lagrangian approach as it exists today occurred 
in 1970 when Held and Karp (1971) used a Lagrangian problem based on 
minimum spanning trees to devise a dramatically successful algorithm for 
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the travelling salesman problem. This work was not only an eye-opening 
successful application, but also set out many key ideas in applying the 

Lagrangian relaxation method to integer programming problems. Fisher 

(1981) and Geoff rion (1974) provide insightful surveys of Lagrangian 
relaxation and its uses in integer programming. Orlin and Bertsimas 

(1991) have developed the most efficient algorithms (in the worst case) for 

many classes of Lagrangian relaxation problems. Belling-Seib, Meyer and 
Muller (1988) have compared Lagrangian relaxation technique for solving 
the constrained network flow problems with other techniques (primal simplex 
algorithm and dual method). 

2.3 LAGRANGIAN RELAXATION TECHNIQUE 

To describe the general form of the Lagrangian relaxation procedure, 
we begin with the following combinatorial optimization model formulated in 

terms of a vector X of a decision variables : 

* 

Z = minimize CX ( p ) 

subject to : 

AX 3 b, (2.2a) 

X € X. (2.2b) 

The model has a linear objective function CX and a set AX ^ b of 
linear explicit linear constraints. The decision variables X are also 
constrained to lie in a given constraint set X = {X : INX = q, 0 ^ x ^ u> 
might be all the feasible solutions to a network flow problem with a 
supply/demand vector q. We further assume that X is finite. 

Lagrangian relaxation procedure uses the idea of relaxing explicit 
constraint set (2.2a) by bringing them into the objective function by 
associating a Lagrangian multiplier vector p. The resulting problem will 
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be: 

minimize (CX + p (AX - b) (L) 

subject to X e X. 

We can refer to the function 

L(#x) = min {CX + p ( AX - b) : X e X>, (2.2c) 

as the Lagrangian function. 

If the constraint set (2.2a) is less than or equal to type, then p 
will be restricted to the non-negative, i.e., ip - 0) values; else, if it 
is equality constraint set (AX = b), then p will be unrestricted in sign; 
otherwise, if constraint set (2.2a) is greater-than or equal to type (AX 
b), then p will be restricted to negative ip < 0). 

Theorem 2.1 : (Lagrangian Bounding Principle). For non-negative Lagrangian 
multipliers (p £ 0), associated with constraint set (2.2a), i.e., AX s b, 

the value L(p) of the Lagrangian function (2.2c) is a lower bound on the 
optimal objective function value Z* of the original optimization problem 
(P). 

Proof : Since AX s b for every feasible solution to (P) and for 

non-negative Lagrangian vector set (|i a 0), the result will be : 

L(fi) * CX* + pi AX* - b) £ (CX* = Z*). 

* * 

We assume that X is an optimal solution to (P). As (AX -b) will be 

* 

either zero or less than zero for an optimal solution X to (P), then for p 
fc 0 , p (AX*-b) will be either zero or negative. Hence Lip) will be a 
lower bound on the optimal value 2 for the original optimization problem. 
To get the sharpest possible bound, we would need to solve the following 


optimization problem : 
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* 

L = max L(p) , 
tiz 0 

which we refer to Lagrangian multiplier problem associated with the 
original optimization problem (P) . The Lagrangian bounding principle has 
the following implication. 

Theorem 2.2: (Weak Duality Theorem). The optimal objective function L of 
the Lagrange multiplier problem is always a lower bound on the optimal 
objective function value of the problem (?), i.e., L* < Z * . m 

Theorem 2.3: If for some choice of the Lagrangian multiplier (p z. 0), the 
solution X of the Lagrangian relaxation (i) is feasible in the 
optimization problem (P), and (ii) satisfies the complementary slackness 
condition p(AX -b) = 0. Then X is an optimal solution of the original 
optimization problem (P). 

Proof: By assumption, L(/i) = CX + p(AX - b). Since p(AX - b) = 0, 

* * 

Lip) = CX and the solution X is feasible to the optimization problem (P), 

* * 

i.e., AX ^ b. Hence, X is an optimal solution to the optimization 

problem (P). m 

The above mentioned properties and theorems show that certain 
solutions of the Lagrangian subproblems provably solve the original problem 
(P). From this, we arrive at two properties : 

Property 2.1: Solutions to the Lagrangian subproblem (L) that are 
feasible, but are not provably optimal for the original problem (P). 

Property 2.2: Solution to the Lagrangian subproblem (L) are not feasible 
to original problem (P). 

In the first case (Property 2.1), branch and bound procedures to get 
the sharpest possible bound can be used. Fisher (1981) provided an insight 
to the use of Lagrangian problem in place of linear programming relaxation 
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to provide bounds in branch and bound algorithms. In the second case, for 
many applications, researchers have been able to modify "modestly'’ 
infeasible solutions so that they become feasible with only a slight 
degradation in the objective function value by using some heuristic 
methods. 

Subgradient Optimization: 

The following issue is related while solving the Lagrangian multiplier 

* 

problem (L): How to find out the optimal multiplier value (p ) of the 
Lagrangian multiplier problem (L), i.e., we need to find the highest point 
of the Lagrangian multiplier function L(p). 

For our optimization model (P) defined as {CX : AX ^ b, X e X>, we 
assume that the set X = {X 1 ,X 2 , . . . 9 )P} is finite. By relaxing the 

constraints AX £ b, we obtain the Lagrangian multiplier function L(p) = min 
{CX + p(AX-b) : XeX}. Thus : 

L(#i) * CX k + p( AX* - b) for all k = 

The best choice for p would be an optimal solution to the dual problem : 

max w (D) 

subject to : 

w ^ CX^ + p(AX^-b), for all k = 1,...,K. 
p i 0. 

The above linear program (D) can be solved by applying linear 
programming methodology. Dantzig-Wolfe decomposition or generalized linear 
programming, is an important solution strategy to solve the linear program 
(D). It has been discussed by Ahuja, Magnanti and Orilin (1993). However, 
one of the disadvantages of this approach that it requires the solution of 
a series of linear programs which are rather expensive computationally. 
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Another approach to find the optimal value of Lagrangian multipliers 

might be to apply gradient method to the Lagrangian function L(/i). The 

function L(p) has all properties like continuity and concavity except one- 

differentiability. The function [L(/i)] is non-differentiable at any fi, 

where it has two or more solutions. For solving the (non-differentiable) 

Lagrangian, subgradient optimization technique is commonly used. The 

subgradient method is a brazen application of the gradient method in which 

gradients are replaced by subgradients. Let /i° be any initial choice of 

k 

the Lagrange multiplier; the subsequent values p for k = 1,2 of the 

Lagrange multipliers can be determined as follows : 

, k+ 1 k „ , . ,Jc , , , + 

[ p = p + S^CAX - b)] . 

k 

In the above expression, X is any solution to Lagrangian subproblem (L) 
k 

when p = p and 6^. is the step length at the kth iteration. As the updated 
value of one or more components of p may become negative, we avoid this 
possibility by taking only the positive value. The notation [Y] denotes 
the "positive part" of the vector y; that is, the ith component of [Y] 
equals the maximum of zero and Y^ . 

Held, Wolfe and Crowder (1974) have discussed the computational 
performance and theoretical convergence properties of subgradient 
optimization in detail. The step size 0 k used commonly in practice is : 

0 = tUB - L(p k )] 

k ||(AX k -b)|| 2 

In this expression, UB is an upper bound on the optimal objective function 
value Z* of the problem (P) and it can be any known feasible solution to 
the problem CP) . 
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Relationship to Linear Programming 

Two properties are important in evaluating a relaxation : the 
sharpness of the bound produced and the amount of computation required to 
obtain these bounds. The primary use of the Lagrangian relaxation 
technique is to obtain lower bounds on the objective function values 
(discrete) optimization problems. By relaxing the integrality constraints 
in the original problem (P), we obtain an alternative method for generating 
lower bound, known as linear relaxation. 

Theorem 2.3: Integrality property : Lagrangian function L(p) satisfies 
the integrality property if it has an integer optimal solution for every 
choice of p even if we relax the integrality restriction on the variable X. 

Fisher (1981) has proved that for problems satisfying the integrality 
property, solving the Lagrangian multiplier problem (L) is equivalent to 
solve the linear relaxation of the problem (P). In these situations, 
though Lagrangian relaxation technique provides no better a bound than the 
linear programming relaxation, yet, the Lagrangian relaxation technique 
might be of considerable value, because solving the Lagrangian multiplier 
problem (L) might be more efficient than solving the linear programming 
relaxation directly. It will be an efficient solution strategy where 
Lagrangian subproblems are network flow problems in nature. The Lagrangian 
relaxation can exploit the core network substructures - shortest path, 
minimum cost flow, spanning tree, assignment and other problems. Usually, 
we choose the constraints to relax so that we can exploit underlying 
structure of the relaxed problem. 

Our work is concentrated to identify all such applications in a 
variety of areas where Lagrangian relaxation of complicating constraints 
results in network flow subproblem (or subproblems). These Lagrangian 
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subproblems can be solved efficiently by network flow algorithms for a 
fixed value of Lagrangian multiplier vector and we can get bounds for the 
optimal solution of the original problem. Ahuja, Magnanti and Orlin (1989) 
presented a literature survey of network flow problems and described some 
of the fastest algorithm for solving these network flow problems. 
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CHAPTER 3 

APPLICATIONS OF SHORTEST PATH PROBLEMS 

In this chapter, we report applications in which the relaxed 
Lagrangian subproblems are shortest path problems. The following is the 
list of applications reported: 

1. Optimal routing of hazardous material vehicle 
(Ramgopalan, Batta and Karwan [1990]) 

2. Routing of a by-product pipe line 

3. Optimal election campaign tour 
(Ramgopalan, Batta and Karwan [1990]) 

4. Optimal routing of a irrigation canal 

5. Routing of a police car trip 
(Ramgopalan, Batta and Karwan [1990]) 

6. Approximating piecewise linear function 
(Imai and Iri [1986]) 

7. Multi-period assignment of consultants 
(Aronson [1986]) 

8. Multi-period assignment of machine in fixed layout problems 

9. Optimal assignment of sales force 

10. Traffic assignment in communication satellites 
(Balas [1983]) 

11. A shortest path problem with side constraints 
(Aneja et al. [1983] and Minoux [1989]) 
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12. School bus scheduling 

(Desrosiers, Sauve and Soumis [1988]) 

Application 3.1 : Optimal routing of hazardous material vehicle 
Reference : Ramgopalan, Batta and Karwan (1990) 

A city municipal corporation has to plan the route of a hazardous 
material vehicle to transport material from a waste collection centre to 
isolated disposal centre. As there is a risk of exposure to city 
population from hazardous material, the corporation wants to minimize the 
risk. Furthermore, city is divided into IK mutually disjoint zones or areas 
labeled as z , ...,z . The objective is to minimize the total risk of 

1 trv 

exposure to the population of the city with an equitable distribution of 
risk among all the zones of the city. 

We define a city transportation network G = (N,A) with N as a set of 

nodes representing the junction of two roads, and A as an arc set 

representing roads of the city. Origin (0) is the waste collection centre 

and destination (D) is the waste disposal site. Between any two nodes i 

and j, arc (i,j) exists if a road connects the junction i and junction j. 

Furthermore, a road may pass through several zones (z^-.-.z^) of the city. 

For each link (i.j) e A and k € C.j * 0 represents the global 

risk associated with travel on arc (i,j), and n (i.j) - 0 represents the 

z k 

risk to zone z, associated with travel on arc (i.j) such that . 

iv 


K 

V - 1 



The corporation wants to have an equitable distribution of the risk 
among the various zones, that is, risk between any two arbitrary zones 
should be under a threshold of p (predeterminded) for the minimum risk path 
of the hazardous material vehicle. 
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The problem can be formulated as an equity constrained shortest path 
problem as follows: 

minimize V Y C. . X 

i j 

subject to: 

E E (rc, (i. j) " (i, j) ) X. . £ n for all a,b e 1, . . . ,K, 

i j a z b 1J 

(3.1a) 

' 1 if i = 0 
0 if i * 0, i * D, 

-1 if i = D, 

(3.1b) 

1 if arc (i,j) is in optimal solution 
X = j (3.1c) 

J 0 otherwise. 

V 

Constraint set (3.1a) ensures equitable distribution of risk between any 

2 

two arbitrary zones. For given K. zones, there will be (K -K) constraints 
of type (3.1a). Constraint set (3.1b) and (3.1c) are shortest path 
constraints. 

The problem is a constrained shortest path problem because of the 

presence of additional constraint set (3.1a). The model can be solved as a 

shortest path problem, if we relax constraint set (3.1a) by assigning 

Lagrange multipliers U (£ 0) to each constraint in this set and forming 

z z, 
a b 

a Lagrangian subproblem as follows: 


{j: (i, j)eA> 


X. . 
ij 


< j: ( j, i)eA} 


Ji 


min L(x,u) = ZD C X + £ £ U (£ E (n (i, j)-nr (i,j))X. -p) 

ij JJ zz, abij a b J 

J a b 

subject to the constraint sets (3.1b) and (3.1c). 

Alternatively, we can say: 

min Ltx.u) = £ E k, ♦ E I U z <” 2 U.j) - H (i,j))] x - Ml E U 2 

i i L J z z ab a b zz.ab 

J a b a b 
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subject to the constraint sets (3.1b) and (3.1 C ) 

For a given vector of multipliers u, £ u ) will be a 

z z. Z a Z b 
a b 

constant. Hence, Lagrangian subproblem can be solved easily by finding a 
shortest path with modified link risks as : 


C'. . 
ij 


(C ij + ^ 1 U z z J n z (i ’j } “ * z (i.J))). 

Z a Z b a b a z b 


such that : 


C( . 0 for feasible solutions, 
ij 

Application 3.2 : Routing of by-product pipeline 

In an Industrial city, a chemical company is engaged in manufacturing 
of products whose production results in certain by-product gases at 
different stages of manufacturing cycle. The company wants to utilize 
these by-products, by transporting these gases to some other plants in the 
city. As these by-products gases can be used as utilities at other plants, 
the company may get some revenue. The company has to get permission from 
state pollution clearance board to lay underground pipeline network to 
connect to the demand plant. 

The company wants to design a underground pipeline network in a way, 
to minimize the total risk of exposure of the gases to the city population, 
as well as an equitable distribution of the risk among all the zones/areas 
of the city. The problem can be formulated as an equity constrained 
shortest path problem like Application 3.1. 

Sometimes, a total investment budget constraint for the construction 
of such a pipeline can also be added. Let a^ be the expected investment 
to construct arc (i,j) and d be the total budget limit for construction. 


Then : 
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I £ a. . X. . < d. 
• j 1J ij 


(3. Id) 


The problem can be solved as a shortest path problem by relaxing total 

investment constraint (3. Id) and equitable distribution constraint (3.1a), 

using Lagrangian multipliers X a 0 and U a 0, with modified arc costs 

z z , 
a b 

as : 


C ij = (c ij + E E U z z. ( ” 2 (1 -» - (1 ' J)) 

^ ^ z z, aba b 

a b 


where C[ . 0 for feasible solutions, 

ij 

This application can be found easily in the 
chemical industries. 


+ Aa. .) . 
ij 


case of process and 


Application 3.3 : Optimal election campaign tour 
Reference : Ramgopalan, Batta and Karwan (1990) 


A political party has to plan the tour of its main leader before the 
mid-term elections. The party wants to maximize his voting potential 
throughout the country. Furthermore, it also wants to maintain an equity 
in the time, that he spends among all the geographical/political zones of 
the country. The country has been divided in several geographical/ 
political zones. 

The optimal election campaign tour problem is having a similar 

structure that of Application 3.1. Each arc length represents the 

gained voting potential while traversing the arc (i,j) of the 

transportation network and n ^ (i,j) is the gained voting potential in zone 

a 

z , if arc (i,j) passes through the zone z , such that : 
a a 

K 
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The objective function will be a maximization in nature, hence the 
resulting problem will be a longest path. 

Application 3.4 : Routing of a irrigation canal 

A state irrigation board of a developing country has to plan the route 

of a proposed irrigation canal to facilitate the irrigation in 

underdeveloped region of the state. The board wants to maximize coverage 

of the agricultural towns while maintaining equitable coverage of the 

districts in that region. The problem can be formulated as Application 

3.1, with the only difference that the objective function will be 

maximization in nature. To avoid looping paths (cycles) in the problem, we 

can insist that on X. . = 0 or 1 , and £ X^ ^ 1. 

ij j J 

Some times, a budget constraint (3. Id) can also be added on the total 
investment limit for the construction of the canal. 

Application 3.5 : Optimal routing of a police car trip 
Reference : Ramgopalan, Batta and Karwan (1990) 

A city police department has to plan the route of a police patrolling 
car trip, to maximize the coverage received by the citizens while equitably 
distributing its time among the communities for which it is responsible. 
The problem is also similar to Application 3.1. 

Application 3.6 : Approximating piecewise linear function 
Reference : Imai and * ri (*986) 

Numerous applications encountered within many different scientific 
fields use piecewise linear functions. These functions contain a large 
number of breakpoints, hence they are expensive to store and manipulate. 
In these situations, it might be advantageous to replace the piecewise 
linear function by another approximating function that uses fewer break 
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points. By approximating the function, we will generally be able to save 
the storage space but we will, however, incur a cost because of the 
accuracy of the approximating function. Assume that, we want to 
approximate the function using a given number p of points. The objective of 
making approximation is, that we would like to make the best possible 
tradeoff between these conflicting costs and benefits such that the new 
function passes through a given number of points. 


Let f^(X) be a piecewise linear function of a scalar X. We represent 

the function in the two dimensional plane; it passes through n points a^ = 

(X Y ), a„ = (X_ , Y 0 ) , . . . ,a = (X ,Y ). Suppose that we have ordered the 
1 1 2 22 n nn 

points so that - x 2 ^ - X^. We assume that the function varies 

linearly between every two consecutive points X i and X i+1 - We consider the 
situation in which n is very large and for practical reasons we wish to 
approximate the function f^(X) by another function f 2 (X) that passes 
through not more than p points {inclduing a 1 and a n >, such that p < n. 

The approximated function results in the savings in storage space. 
Let a be the per unit cost associated with any single interval used for 
approximation. Furthermore, the error of an approximation is proportional 
to the sum of the squared errors between the actual data points and the 

r J ; 

estimated data points, i.e., the penalty is 0 £ (f^X^ - f 2 ^ X k^ 


for 


some constant /3 . The problem can be formulated as a constrained shortest 
path problem on a network G = (N,A) with n nodes, numbered 1 through n. 
The arc (i,j) signifies that we approximate the linear segment of the 
function fj(X) between the points a., a. +1 ,...,a j by one linear segment 
joining the points a. and a y The cost of the arc (i, j) will be given by : 

c jj - « ♦ p [J. <W - W )2 ]' 




Fig 2 : Illustrating Application 3.6 : (a) approximating the function fjCX) 
passing through 10 points by the function f 2 (X) ’ ^ corresponding 
shortest path problem. 
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Let X iJ be a zero-one variable, indicating whether arc (i.j) is included or 
not in the optimal solution. The problem can be formulated as : 


minimize T C X 
(1. j)sA 1J 


subject to: 


E x n ' E X.. 

{J:(i,j}eA> J {j:(j,i)eA} J1 


E X s p-1 for all i, j e N, 
(i,j)€A 1J 


+1 for i = 1, 

0 for i * l,n, 

-1 for i = n, 

(3.6a) 


(3.6b) 


( 1, if arc (i.j) is included, 
0, otherwise. 


(3.6c) 


Constraint (3.6b) ensures that the function should be approximated in 
a manner, such that there should not be more than p points in the new 
function. This constraint can be relaxed, using a Lagrangian multiplier A 
i 0. The relaxed problem will be 


minimize L(x,A) = T (C. . + A) X. . - A(p-l). 

(i,J)€A 1J 

subject to the constraint set (3.6a) and (3.6c). 


For a fixed value of A, A(p-l) will be a constant. The relaxed 

problem can be solved as a shortest path problem having modified arc costs 

as C' = C. . + A. Note that here, we assume that the new function must 
ij iJ 

pass through the points 1 and n. 

Application 3.7 : Multiperiod assignment of consultants 
Reference : Aronson (1986) 

A consultancy firm provides project monitoring and project assistance 
services to its clients. The firm appoints its consultants at different 
client offices depending upon the project requirements. The firm wants to 
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determine the optimal assignment of its consultants over a fixed time 

period T. With each assignment of consultant i to client j in any time 

period t (t = the assignment cost. The assignment cost 

need not to be equal in each time period because of varying skills of 

consultants and different requirements of projects. Furthermore, let S , 

i jkt 

be the transferring cost (training, moving, resettling) of consultant i 
from client j in period t to client k in period t+1. 

The optimal multiperiod assignment of consultants problem can be 
described as an integer multicommodity flow problem. We define a network G 
= (N,A), consisting of a set of nodes N and a set A of ordered pair of 
nodes (p,s), called arcs. Each consultant is represented as a commodity. 
We also assume that the number of consultants and number of clients are 
both equal to n. 

For the network G = (N,A), let N + £ N denote the set of nodes s 6 N 

P 

for which arc set (p,s) € A, i.e., the set of nodes which have arcs 

pointing away from node p; and Q N denote the set of nodes s € N for 

which the arc (s,p) e A, i.e, the set of nodes which have arcs pointing 

toward node p. Furthermore, for each commodity i, i = l,...,n, with each 

arc ( p , s ) € A is associated a cost per unit flow C* which will be set to 

P s 

either C. or S. . . We define the node set to be N 

ljt ljkt 

{0,1, . . . , (2T-2)n+l)> and arc set to be the union of 2T-1 disjoint sets; T 
sets of assignment arcs (A^.) and T-l sets of transfer arcs (T^). The node 

requirements for i = l,...,n are r* = +1, r (2T-2)n+l = anc * r p = 0 for 

p = 1, . . . , (2T-2)n. The problem can be formulated as follows: 

n i i 

minimize £ £ C X 

i=l (p, s)eA ps PS 
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subject to 




£ - X sp ~ V P e N, i 
seN sp p 
P 


1. • • • ,n. 


n 


E 

i=l 



X ps “ (p ’ s) € A t - t = 1 T, 

i 0, and integer, (p,s) e A, i = 1 n . 


(3.7a) 

(3.7b) 

(3.7c) 


The mutual capacity constraints (3.7b) prohibits the assignment of a 
client to more than one consultant in each period. We propose to relax the 
mutual capacity constraints (3.7b) by using nT Lagrangian multipliers 
(Ap^ £ 0). The relaxed problem will be to 


i i n 

minimize £ £ C X + £ A ( £ X 1 - 1) 

i=l (p,f)€A PS PS ( Pt s)eA t PS 1=1 ps 


subject to (3.7a) and (3.7c). 


or* equivalently, 


n 


n 

min £ £ (C 1 + A ) X 1 + £ £ ^ A - , 

1=1 (p.s) 6 A t ps ps ps lil (p,s)sT t ps ps (p.s)*A t ps 


c 1 x 1 


subject to (3.7a) and (3.7c). 


The resulting objective function for the Lagrangian subproblem has a 

i i 

cost of C + A for each assignment arc and C for each transfer arc. 

ps ps ps 

Since none of the constraints in this problem contains the flow variables 

(X 1 ) for more than one of the consultants, the problem can be decomposed 

ps 

into n pure independent shortest path subproblems. 

Application 3.8 : Multi-period assignment of machines in 

fixed-layout problems 

A manufacturing division has n machines, having different efficiencies 
or operating costs (C^ of Application 3.7) for different sets of task 
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requirements at n different locations over a time period of T. The tasks at 
various locations require similar type of work but with varying efficiency 
requirements and other work parameters. Furthermore, any machine can be 
operated at any location to perform the required task but with varying 
efficiency and costs. In addition, the cost of moving or transferring 
^i jkt ^PP^* ca tion 3.7) these heavy machines from one task location to 
another location may be substantial. The objective of the problem is to 
find out a minimum cost assignment of the machines at different tasks 
locations over a time period T. The problem of optimal multi-period 
assignment of machines in fixed layout problems is also similar to 
Application 3. / . This problem can also be found easily at construction 

sites where mobile machinery is used at different locations and at docks 

where ship construction work is going on. 

Application 3.9 : Optimal assignment of sales force 

A marketing division of an organization wants to determine the 
minimum cost assignment of sales-managers to look after sales and customer 
services at different marketing territories over a time period of T. Based 
on good demand forecasting methods, they know job requirements and work 
load in these territories during each time period in advance. Furthermore, 
sales manager have varying marketing skills and different assignment costs. 
The optimal assignment of sales force problem aims at determining sales 

managers assignment to minimize the total assignment and transfer costs 

over a fixed period of time. 

The problem is having a structure similar to that of Application 3.7. 

Application 3.10 : Traffic assignment in communication satellites 
Reference : Balas (1983) 

Everyday a rapidly increasing volume of long distance TV, radio and 
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telephone communications is transmitted digitally via satellites. A high 
capacity communication satellite interconnects simultaneously scores of 
transmitting and receiving stations. One of the advanced techniques for 
operating such a satellite communication network is SS/TDMA 
(satellite-switched time division multiple access system) based on the use 
of highly directed spot beam antennas. Under this scheme, each transponder 
on board of satellite is allocated to a pair of ground stations for certain 
amount of time, based on the communication traffic between these two 
stations. Each transponder has different capacity. Such a set of 
allocations involving all the transponders is called a 'switch mode'. 
These allocations that make up a switch mode are changed simultaneously by 
an on-board switching facility. The whole sequence of switch modes that 
make up the schedule of allocation for a given period of time is called a 
frame. 

The problem of optimal traffic assignment of communication satellite 
deals with the efficient scheduling of time slot allocations for a frame. 
As each transponder is having different capacity, the problem is to find a 
multi-period assignment of all transponders to various pairs of ground 
stations, based on the communication traffic. 

The problem can be formulated similar to Application 3.7. 

Application 3.11 : A shortest path problem with side constraints 

In many real life problems we have to identify a shortest path from 
one node to another node with an additional constraint. This problem can 
be viewed as multiple criteria shortest path problem. Below, is the list 

of such applications : 
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C A ) A minimum travel cost path between two cities subject to 
constraints on traversal time, energy consumption, 

[Aneja et al. (1983)]. 

(B) Shortest path with fading out constraints for routing problems in 
telephone networks [Minoux (1989)]. 

(C) Path of average minimal length with constraints on the variance 
in probabilistic graphs [Minoux (1989)]. 

These problems can be formulated as shortest path problems having a 
side constraint on the other criteria, as follows : 


minimize Y C. . X 

(i, j)£A iJ 1J 


subject to 


l 

j€N 



I X 

j€N 


ji 


1 if i = source, 

0 if i * source, i * sink, 
-1 if i = sink 


i 

(i, j)6A 




s b, 


X, , € (0,1). 
i J 


(3.11a) 

(3.11b) 

(3.11c) 


The constraint. (3.11b) is a side constraint in the shortest path 
problem. The problem can be solved as a shortest path problem, by relaxing 
complicating constraint (3.11b), using a Lagrangian multiplier X > 0. The 

Lagi angian subproblem will be as follows : 


minimize 


(C ; ; + X a. .) X - Xb 


1J 


ij ij 


(i J)eA 

subject to the constraints (3.11a), (3.11c). For a fixed value of 

multiplier A * 0, Ab will be a constant. The relaxed problem can be solved 


as a shortest path problem. 
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Application 3.12 : School bus scheduling problem 
Reference : Desrosiers, Sauve and Soumis (1988) 

School buses have to be scheduled over a set of trips. In the 
morning, a bus picks up students along its trip and delivers them to the 
school. In the afternoon, a trip begins at the school and ends at the stop 
where the last student is dropped off. Furthermore, it is also desirable 
to have a single bus sequentially cover several trips. Each locality where 
the students reside and every school has to be visited within a specified 
time interval. The problem of school bus scheduling aims at finding the 
minimum number of buses required to visit a set of nodes subject to the 
time window constraint, for a homogeneous fleet of buses located at a 
common depot. This problem was formulated as a network flow problem with 
additional time window constraints by Dersosiers et al. (1988). 

Two I.agiangi.m relaxations of the model were studied by the Dersosiers 
(1'iKH). In the first one, the time constraints are relaxed producing 
shortest path subproblems which are easy to solve, but the bound obtained 
is weak. In the second relaxation, constraints requiring each node to be 
visited are relaxed producing shortest path subproblems with time window 
constraints and integrality conditions. 
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CHAPTER 4 

APPLICATIONS OF MINIMUM SPANNING TREE PROBLEMS 

In this chapter we report applications in which the relaxed Lagrangian 
subproblems are either a undirected minimum spanning tree problem or a 
directed minimum spanning tree problem. We describe the following 
a PPlications in this chapter: 

1. Optimal drilling of oil wells 
(Akinc and Srikanth [1992]} 

2. Optimal tour of a health care delivery team 
(Akinc and Srikanth [1992]) 

3. Optimal tour of a maintenance team 
(Akinc and Srikanth [1992]) 

4. Tour of a mobile departmental store 

5. Travelling buyer’s problem 
(Akinc and Srikanth [1992]) 

6. Servicing of a wide forest area 
(Akinc and Srikanth [1992]) 

7. Design of centralized computer network 
(Volgenant [1989] and Gavish [1982]) 

8* Multi-drop terminal layout in teleprocessing design 
(Boorstyn and Frank [1985]) 

9. Reconstructing of a utility network after natural 
disaster 
(Shogan [1983]) 
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10. Optimal offshore pipeline network design 

(Rathforb, Frank, Rosenbaum and Steglits [1970]) 

Application 4.1 : Optimal drilling of oil wells 
Reference : Akinc and Srikanth (1992) 

Oil and gas commission of a country wants to drill exploratory oil 
wells in a newly discovered oil field to determine the size and other 
properties of oil field. To drill such exploratory wells, platforms are 
needed at these locations for mounting of oil rigs. Furthermore, it is 
usually possible to drill several such wells from a platform that is 
stationary at a point. As the construction cost of these platforms are 
quite substantial; therefore, significant economies may be realized by 
using towable platorms rather than stationary ones and drilling of 
oil-wells from these platforms. The drilling cost will be related to the 
distance and the angle of well to be drilled from the platform. Thus a 
towable platform will act like a mobile service unit that will visit some 
oil well locations physically, but not all the locations, and we will drill 
the rest of the locations from afar (at a cost) while stationed at a nearby 
locations. The OGC wants to determine the subset of all n oil well 
locations the towable platform should actually visit, the sequence of these 
visits and subset of all the points the platform should drill from each 
location on its path so that the total cost of drilling is minimized. 

Let M denote the cost of moving of a towable platform from location 

ij 

i to location j (sum of towing, set-up and take down cost of the platform) 

and P. . denote the cost of drilling at location j while towable platform is 
ij 

at location i. Furthermore, m. denotes the set of oil well locations the 
platform can move to from location i and p i denotes the set of oil well 
locations that can be drilled from location i. S is the set of all the n 
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locations of oil wells to be drilled. 


Let us define a zero-one variable X_ indicating whether or not the 
platform moves from location i to location j, and Y_ indicating whether a 
location j is serviced by the platform stationed at location i or not. We 
add one dummy location 0 in the set S, such that 


S' s S u 0. 

Furthermore, M_ . = M. = 0 ,for all i € S and P^= P. n = for all 
Oi 1 0 Oi lO 

i 6 S. 


The problem can now be formulated as : 


minimize £ Z 
i j€S' 


M. . X. . 
ij 1J 


+ e i 

i jeS' 


P. . 
ij 


Y. . 
IJ 


subject to 


Z Y - |p | E X so 

jeS 1J JeS' 1J 




i*j 


ieS' 


(Y. . + X. .) = 1 
ij iJ 


E X - E X- - = o 


ieS 


JeS' 


ij 


E x rt . = E x = l , 
As 0j jis j0 


Z E (X + 

- J 


j€S k i€S k 


• Y. . s IS. 

J ij 1 k 1 


X. Y. . = 0, 1, 
ij iJ 


for all i e S, 

for all j e S, 

for all i e S, 

1 for S. £ S, 
k 

for (i, j) € S' , 


(4. la) 

(4.1b) 

(4.1c) 

(4. Id) 

(4. le) 

(4. If) 


Constraint set (4.1a) prevents any drilling from a location unless the 
towable platform visits it. Constraint set (4.1b) requires that each 
location either be visited by the towable platform or else drilled from 
another location. Constraint sets (4.1c), (4. Id) and (4.1e) together 
ensure that the path of the visits by towable platform starts and finishes 
at dummy location 0, is contiguous, and is free of any subtours. 
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The problem can be solved by relaxing the constraint set (4.1a) and 
(4.1c), using two Lagrangian multipliers 0 and p. unrestricted in 
sign, as follows : 


L(A,p) = min £ £ M X + £ £ P Y + £ X (|p| £ X - £ Y 

i icC' • L J i j^ e 5 1 1 IrC' 1 


.) 


i jeS 


i jeS 


jeS' jeS 


ij 






I Mi f I X - I X ) 
ieS \jeS' J jeS' J ' 


ieS v jeS 

subject to the constraints (4.1b), (4. Id), (4.1e), and (4. If) 

Alternatively, we can say : 


min I Z 
i jeS' 


L 1J J '‘J 1 


m. . + x. ip. | + n. - p. x. . + y y (p. . - x.) y. . 
’ 1 ’■ i jj ij r 1 


i jeS 


subject to the constraints (4.1b), (4. Id), (4.1e), and (4. If). 

We can also write the objective function as follows 
L(X,p) = min Y V M. . X. . + Y Y P. . Y. . 

ijis' 1J 1J tjts 1J 1J 

where M. . = (M. . + X. Ip. I + u. - u.) and 
ij iJ J j i J 

P. . = (P. . - X. ). 

IJ iJ i 

The optimal solution of this problem has the following properties : 

(a) for any (i,j), X„ and Y„ can not be equal to one simultaneously due 
to the constraint set (4.1b); 


(b) P. . > M. . ==> Y. . = 0; and 

ij ij iJ 

(c) M. . > P. . ==> X. . = 0. 

ij ij ij 


The Lagrangian subproblem can therefore be simplified further by 


defining a new set of variable Z 


where 
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with associated objective function coefficients as : 


C. . = min {P. M. 
ij ij 1J 


For a given set of and the simplified Lagrangian subproblem will be: 

Lin.n) = min £ £ C Z 

i jeS' J J 
i*j 

subject to 


£ z. . = 
iJs' 1J 

1 

for j e S 

(4. lg) 

i*j 




l Z 0j = 

jeS' J 

1 Z j0 

jeS' JU 

= 1 

(4. lh) 

I 2- • 
i.J*s k 
i*j 

£ l S kl - 

1 for S. £ S 

k 

(4. li) 

Z. . = 0 or 

1 

for all i,j. 

(4. Ij) 


The optimal solution to the Lagrangian subproblem (L) will be a set of 
n+1 arcs consisting of the n arcs in the directed minimal spanning tree on 
the nodes in S, with the root node 0 restricted to degree 1, plus a single 
arc returning to node 0. The problem can be solved in polynomial time. 
(Tarjan [1977]; Lawler [1976]). 


Application 4.2 : Optimal tour of a health care delivery team 
Reference : Akinc and Srikanth (1992) 

Health department of a underdeveloped country has to plan the route of 
a mobile health care delivery team in a particular region of the state. 
The team will make a tour of subset of all the villages in the region and 
those needing care at other villages may come for service while the team is 
temporarily stationed at some close-by village. The problem can be modeled 
as in Application 4.1. 
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Application 4.3 : Optimal tour of a maintenance team 
Reference : Akinc and Srikanth (1992) 

A maintenance team has to render routine service of maintaining a 
communication system at geographically dispersed points. As the total cost 
of visiting each point increases, significant economies may be obtained by 
organizing a small maintenance unit (or units) and co-ordinating some of 
their resource needs at temporary and movable bases stationed at some 
nearby locations. These bases may also schedule the activities of small 
units, house hard-to-move equipments, maintain communication with 
headquarters and provide lodging. 

The maintenance team wants to determine the subset of all location 
that it should physically visit, the sequence in which these visits should 
be made and subset of all locations the maintenance team should service 
from each location on its path so that the total cost is minimized. 

The application is having structure similar to that in Application 4.1. 
Application 4.4 : Tour of a mobile departmental store 

A consumer co-operative society has to plan the route of a mobile 
departmental store through various parts of the city. Because of the time 
constraint, the store can’t visit all the localities physically rather it 
will visit a subset of all localities of the city and customers needing 
some consumer goods at other localities may come for purchasing while store 
is stationed temporarily at some close-by locality. The problem is also 
similar to Application 4.1. The cost of providing customer services ( 
of Application 4.1) indirectly from some other nearby locality can be 
measured as expected sales loss because of not serving these customers 


directly. 
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Application 4.5 : Travelling buyer's problem 
Reference : Akinc and Srikanth (1992) 

A travelling buyer has to purchase a number of items from various 
vendors located at different geographical points. These vendors stock a 
subset of all the items and sell them at varying prices. The problem is to 
determine the subset of vendors and the sequence of the visits to these by 
the buyer so that the sum of travel and purchase cost is minimized. The 
traveling buyer’s tour problem can be modeled as Application 4.1. The tour 
of the buyer may be modeled as the moves of a mobile service unit (buyer), 
while the items to be bought while at some vendor’s location may be 
considered as services to be performed. 

Application 4.6 : Servicing of a wide forest area 
Reference = Akinc and Srikanth (1992) 

A forest department has to plan the optimal servicing plan of a wide 
forest area. The services are provided as a routine service by a team of 
forestry personnel to the remote parts of the forest. In rendering these 
services, economies may be obtained by a "satellite" service unit (or 
units) and coordinating their needs for equipment etc. from temporary and 
movable bases. The problem is similar to Applications 4.1 and 4.2. 

Application 4.7 : Design of centralized computer network 
Reference : Volgenant (1989) amd Gavish (1982) 

In a centralized computer network, several terminals at various sites 

have to be connected to a central computer, with a minimum amount of wire. 

Furthermore, the size of the terminals is such that the number of wires 

incident to a terminal k can not exceed a number r Let C. . be the amount 

K 1 J 

of wire required to connect two terminals i and j. Furthermore, we define 
a set of zero-one variables X„, which attains the value of one if arc 
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(i, j) is included and zero otherwise; and Y. . is a set of variable that 

iJ 

specifies flow on the arc connecting terminals i and j. 

The problem of designing an optimal centralized computer network will 
be a degree constrained spanning tree with a root in the central site that 
spans over the other sites and with each terminal k restricted to a degree 
r^ and has a minimal connection wire length. Let I = {2,...,n> denote the 
set of different terminal locations and the central computer site is 
denoted as node 1, such that 

I = I u 1. 


The problem can be formulated as : 

n n 

minimize z = £ £ 

i=2 j=l 
j*i 


C. . 
ij 



subject to 


n 

E X = 1 for i = 2 n, 

j=l J 
j*i 

n n 

E Y. . - E Y . . = 1 for I * 2 n, 

J=1 1J j=2 J1 


(4.7a) 


(4.7b) 


Y. . ^ (n-1 ) X. . 
1J ij 


for i = 2, . . . , n, 

j = 1,2 n and i * j, (4.7c) 


Y. . i 0 and X. . = 0 or 1 
ij ij 


n 


n 


l x lk - l x ki s 

i=2 llc i=I K1 K 


for all i,j e I, 

for all k e S. 


(4. 7d) 

(4. 7e) 


Constraint sets (4.7a), (4.7b), (4.7c) and (4.7d) ensure the spanning tree 
structure of the network and constraint (4.7e) ensures that for any 
terminal k, its degree should not exceed r^ i 1. The subset S is a set of 
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all nodes on which degree constraint has been imposed. 


The problem can be solved as a minimal spanning tree problem by 
multiplying the constraints in (4.7e) by a vector of Lagrangian multiplier 
{A} and adding them to the objective function. This gives the following 
Lagrangian subproblem : 


n n 

L(A) = min £ A r + £ £ C X 

keS 1=2 j=l J J 


subject to the constraints (4.7a), (4.7b), (4.7c) and (4.7d), where 


C. . 
1J 


for i g S and j g S, 

C. . - A. for i € S and j £ S, 
ij i J 

C..-A. for i ^ S and j € S, 
iJ J 

C. . - A. -A. for i € S and j € S, 
1 J 1 J J 


For a fixed vector of multiplier A, the problem will be a minimal spanning 
tree problem. 

Application 4.8 : Multi-drop terminal layout in teleprocessing design 
Reference : Boorstyn and Frank (1985) 

Computer communication networks are becoming a prevalent part in the 
corporate environment. A problem, which is common to many computer 
communication network design problem is to find a minimal cost spanning 
tree that connects a set of known terminal sites to the central site. In 
such networks, generally a polling scheme to regulate the traffic is used 
so that the terminals can send messages over common branches. The tree 
structure of a teleprocessing network can be divided into subtrees rooted 
at the central site. Each subtree is called a multidrop line. The 
terminals on each multidrop line are regulated by a polling regimen. 
Furthermore, each port of the communication controller in the central site 
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that deals with a multidrop line ha^ ~ 

ds a limited capacity on the amount of 

traffic that it can handle. So. the ^ ^ , . , 

capacity of the link connecting each 

multi-dropline to the central site mu*.*. u 

^ xte m ^st be less than or equal to the 

limited capacity of each port. The ^ , . , . J t , 

problem was formulated as a capacitated 

minimum spanning tree problem by Gavich noon * » . , 

^ a y <iv:Ls n [1985]. A Lagrangian relaxation 

was used to solve the problem as a deerp^ . . , . . . 

& ree constrained minimum spanning tree 

problem. 


Application 4.9 : Reconstructing of a 
natural disaster 

Reference : Shogan (1983) 


utility network after 


Cities, states, and countries rely heaylly on a yarlety of utlllty 

networks : energy networks (such as electrical or natural-gas networks), 

communication networks (such as . , , 

telephone, telegraph, or computer 

networks), transportation networks (such =0 v- u . . . t , , 

'• i: ' ucn as highway or railrood networks), 

and water networks (such as networks f or the dlstrlbutlon of portable „ ater 

or the removal of sewage). Many important networks do not have a tree 

structure, and, after sustaining damage from some natural disaster (such 

as from earthquake), the long run goal for such a netuork „ ln be to 

rebuild the network so that it is as pcoh w ^ .v - 4 . . 

s good or better than it was prior to 

the earthquake. However, a less ambit Inno , , , . 

H u - LL1 °us short run goal may simply be to 

allocate scarce resources (Labor, equipment, spare parts etc.) to the 

repair process in such a way as to restore some -minimal level of service" 

as quickly and economically as possible. 


Let N be the set of nodes- having indices 1>2 ,..., n; node 1 represent 
the source having an infinite supply of a commodity, and every other node p 
is a sink node having a known demand for that commodity. There are m = \ 
n(n-l) arcs in the network before the natural disaster. Associated with 
each arc j(l s j < m), ej is the flow capacity of arc j. The minimal level 
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of service after natural disaster can be provided by constructing a 

spanning tree (T) when only a single source exists, or, by a spanning 

forest (when multiple sources exist); so that every demand node would have 

access to a source. Furthermore, the restortion of arc j(l 5 j < m) to its 

former capacity e. would cost c. dollars and would consume a. . (1 ^ i ^ K) 
y J J ij 

units of the limited quantity of b^ units of resource i. The problem can 
be formulated as a resource constrained spanning tree problem, as follows : 
m 


min Y C. X. 
j-i J J 


subject to 


e 

j=i 


a. . X. s b. 
ij J i 


X € T. 


1 £ i £ K, 


(4.9a) 

(4.9b) 


Here, X denotes a binary vector whose jth component X^ equals 1 if arc j is 
included in the minimal spanning tree and equals 0 otherwise. The problem 
can be solved as a minimal spanning tree problem by relaxing the constraint 
set (4.9a), and associating non-negative Lagrangian vector X = 
(A , A„ , . . . , A„) with each constraint. The relaxed problem will be : 

1 cL fv 

m Km 


L(A) 


mm 


j=i 


C. X. + 
J J 


E \ ( E 


i=l 


j=l 


a. . X. - b. ) 
IJ J 1 


subject to (4.9b), 


or, equivalently. 


m 

minimize £ 

j=l 


C. + Y A. a. . 


l J 


i=l 




X. 

J 



A. b. 


l l 


subject to (4.9b). 

For a fixed value of multipliers A., 


K 

y A. b. will be a constant and the 
“.ii 


subproblem will be a minimum spanning tree. 
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Application 4.10 : Optimal design of offshore pipeline network 
Reference : Rathforb, Frank, Rosenbaum and Steglits (1970) 


Several offshore natual gas fields have to be connected to an onshore 
separation and compression plant. The objective is to find out a mini^ um 
cost spanning tree structure that will connect all the fields to the 
onshore plant. Sometimes, a constraint on the amount of pipe available can 
be imposed in such a design. The problem can be formulated as a resource 
constrained minimal spanning tree problem like Application 4.9. A detailed 
model having certain other constraints on several other criterias has been 
discussed by B. Rathforb et al. (1970). The relaxed subproblem can be 


solved as a spanning tree problem, as in Application 4.9. 
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CHAPTER 5 

ADDITIONAL APPLICATIONS 

In this chapter, we report applications whose Lagrangian relaxations 
result in more than one Lagrangian subproblems . In addition to these, we 
will also report the applications in which relaxed Lagrangian subproblems 
will be network flow problems other than shortest path and minimum spanning 
tree problems. These applications are as follows: 

1. Optimal assignment of trainees 

2. Core management of nuclear reactor 

(Gupta and Sharma [1981]) 

3. Job shop loading 

4. Tour of a news paper vehicle 

5. Optimal tour of a courier service personnel' 

6. Sequencing of a drilling job 

7. Interconnection of local area networks 
(Fetterolf and Anandlingam [1991]) 

8. Electrical distribution system design 
(Current and Pirkul [1991] ) 

9. Optimal design of mass transit system 
(Current and Pirkul [1991]) 

10. Optimal design of gas field pipeline network 
(Current and Pirkul [1991]) 

11. Transportation network design for resource extraction area 
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Transportation network design for a developing country 
(Current, Re~velle and Cohon [1986]) 

13. Optimal newspaper distribution network 
(Current, Re-velle and Cohon [1986]) 

14. Project scheduling with resource constraints 
(Christof ides, Alvarez-valdes and Tamirit [1987]) 

15. Procurement of aviation fuels 
(Austin and Hogan [1976]) 

16. A forest industry problem 
(Jornsten and Vabrand [1986]) 

17. A logistics planning system problem 
(Klingman, Mote and Philips [1986]) 

Application 5.1 : Optimal assignment of trainees 

A heavy engineering company has to recruit workers as 

trainee-apprentice for one of its manufacturing division. The recruited 

workers have to be provided training for different tasks in that 

manufacturing division. The company has some yearly training budget to be 

spent on the training of new employees. Let C_ be the expected revenue 

for the company and be the training expenditure, if ith trainee is 

assigned to jth job. Furthermore, the total number of trainees, n, are 

equal to the number of tasks to be assigned. The numbers C. . and a. . are 

ij ij 

known in advance because of the qualification and technical skills of the 
trainees. The optimal assignment of trainees problem can be formulated as a 
constrained assignment problem to find out the one-to-one assignment of 
trainees and jobs, with a view to maximize the expected revenue of the 
company. In case, the training budget should not exceed more than d 
rupees, this problem can be modeled as follows: 
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n n 

maximize £ £ 

i=l j=l 



subject to 


n 

E X = 1 for j = 

i=l J 

n 

I X = 1 for i = 1, n, 

j=l J 


n n 

I E 

i=l j=l 


a. . 
ij 


X. . 
ij 


( 1 , if ith worker is assigned on job j, 
0 , otherwise. 


(5. la) 

(5. lb) 

(5.1c) 

(5. Id) 


Constraint set (5.1a) ensures the assignment of only one trainee for 
each task. Constraint set (5.1b) ensures assignment of only one task for 
each trainee. Constraint (5.1c) ensures that training expenditure should 
not exceed more than d rupees for the assignment. 


The problem can be solved as an assignment problem, if we relax 
constraint (5.1c) by associating a Lagrangian multiplier A 2: 0 with it. 
The relaxed problem will be as follows : 
n n 


max E E 
i=l j=l 


(C. . + Aa. .) X. . - Ad. 
IJ IJ IJ 


subject to constraint set (5.1a), (5.1b) and (5. Id). 


Application 5.2 : Core management of nuclear reactor 
Reference : Gupta and Sharma (1981) 

A nuclear power plant that uses pressurized water reactors (PWR) has 
to find out the optimal reactivity allocation of the given fuel assemblies 
having different burn up exposures. Let there be n fuel assemblies to be 

assigned to n possible locations in the reactor i, EfffiRWjedn?«AR>-he 

t 1 r .AVPufi 
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problem is to find out the one to one assignment of fuel assemblies and 
reactor locations to maximize its reactivity, under the constraint imposed 
upon the peak load factor of power distribution. Let be the maximum core 
power distribution factor of the reactor. The problem can be formulated as 
follows : 

n n 

iii A 1J 1J 

subject to 


n 

n 

i 

i 

x=i 

j=i 

n 

i 

i=l 

X. . 
1J 

n 


X. . 

j=i 

1J 

X. . 
1J 

■( 


for j = 1 n, 


for i = 1 , . . . ,n, 


(5.2a) 


(5.2b) 


(5.2c) 


1 if fuel assembly i is assigned to reactor location j, 
0 otherwise. 

(5. 2d) 


Here, CL ^ and d_ are non-negative parameters. These can be calculated as 


follows: 


E. 

l 


= exposure of fuel assembly i, 

= core power distribution form factor, 


( Vi 

= Statistical weight 

c Vi 

= Statistical weight 

5 iJ 

= (E.-Ej) 

(w R )i . 

d ij 

= (E. - EJ 

1 J 

(w F )i , 

C ij 

= max C. . - 

i.j 1J 

c. . , 

1J 

d ij 

* 3 u " 

d. . , 
ij 
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= maximum power distribution form factor 
= F + n (min d_ ) . 

This problem is a constrained assignment problem, due to the 
constraint (5.2a). This can also be solved by using Lagrangian relaxation, 
as in Application 5.1. 

Application 5.3 : Job shop loading 

A job shop has a small group of similar machines having varying 
capacity etc. The functional requirement of various machines within a 
group is same, though their processing time, capacity, power consumption 
vary. Due to the nature of job shop, each such group of machines receives 
workpieces of different sizes and shapes in the lots. Let C be the cost 
of assigning workpiece i to machine j and t_ be the processing time for 
such assignment. The job shop loading problem aims to minimize the cost of 
assignment such that the total processing time of each lot should be within 
available ’ d* machine hours. 

The problem is similar to Application 5.1 except that the objective 
function will be minimization in nature. 

Application 5.4 : Tour of newspaper vehicle 

A newspaper printing plant has to determine the minimum distribution 
cost tour of its newspaper vehicle that will transport and deliver 
newspapers to certain cities of a particular zone, near to the printing 
plant. 

The vehicle will start from its home base (printing plant), node 1, 
visit each of several cities in that zone, represented by nodes 2,...,n, 
exactly once and return home base after delivering newspapers, doing so at 
the lowest possible travel cost. Furthermore, to ensure good delivery, a 
time constraint on the total travel time of the vehicle is also imposed, 
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time constraint on the total travel time of the vehicle is also imposed, 
so that demand nodes can get newspapers as early as possible. 


The problem can be formulated as travelling salesman problem with an 

additional constraint on the total travel time of newspaper vehicle. Let 

C. . denote the cost of traveling and t. . denote the travel time from city i 
ij ij 

to city j. Furthermore, plant wants that vehicle trip should not be more 
than T hours. Let us define a zero-one variable indicating whether or 

not the vehicle travels from city i to city j, a flow variable on each 
arc (i,j) and assume that vehicle has n-1 units available at node 1 (home 
base) and it has to deliver 1 unit to each of the other nodes. The problem 
can be formulated as : 


minimize 


I C. - 
(i.J)€A 1J 



subject to 


i 


l-j-n 


Y. . 
IJ 


1 


I 

l^i^n 


Y. . 
ij 


1 


Y X. . - Y X.. = 

J 1J J J1 


X. . * (n-1) Y. . 
ij ij 


for all i = 1,2, ... ,n, 

for all j = 1 , 2, . . . , n, 

n-1 for i = 1, 

-1 for i * 1 , 

for all (i,j) e A, 


t. . Y, . 
ij iJ 


T, 


X. . & 0 and integer, for all (i,j) e A, 
^ J 


Y 


ij 


1, if vehicle visits arc (i,j), 
0, otherwise 


(5.4a) 

(5.4b) 

(5.4c) 

(5. 4d) 

(5. 4e) 

(5. 4f ) 
(5. 4g) 


The constraint (5.4a) and (5.4b) ensures that exactly one arc leaves 
and enters any node i, i.e., integer solution satisfying (5.4a) and (5.4b) 
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will be the union of disjoint cycles. Constraint (5.4c) is flow balance 
constraint and ensures connectivity of the network, constraint (5.4d) are 
redundant if Y_ = 1 as no arc need ever carry more than (n~l) units of 
flow. Constraint (5.4e) ensures that the total trip time of vehicle should 
not exceed more than T hours. 


The problem can be solved by attaching Lagrangian multipliers p_ £ 0 
with the constraint (5.4d) and A with constraint (5.4e) and bring them into 
objective function, giving the Lagrangian objective function : 

min £ [C. . - (n-1) p. . + At. .] Y. . + Y p. . X. . - AT 

Ci.J)€A 1J 1J 1J 1J CI.J)€A 1J 1J 

subject to (5.4a), (5.4b), (5.4c), (5.4f) and (5.4g). 

The relaxed problem decomposes into two subproblems : 

(i) an assignment subproblem in variable Y_ 
min Y [C. . - (n-1) p. . + A t.J Y. . 
subject to (5.4a), (5.4b) and (5.4g). 

(ii) a minimum cost flow problem 


mm 


p. . X. . 
ij ij 


( i , j ) € A 

subject to (5.4c) and (5.4f), 


For fixed value of A, AT will be constant. 

These subproblems are network flow problems. 

Application 5*5 : Optimal tour of courier service personnel 

A courier service personnel has to deliver various packets of letters 
at several cities, collected from the home city. The courier service wants 
to find out minimum cost travel tour, with a time constraint on the 
travelling done by employees. This problem will be similar to Application 
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4.1. Here, the cities, to be visited by messengers are connected through 
transportation network and are in the same zone so that the same person can 
travel through all the cities easily. Such problem may arise when different 
transportation modes are available between two cities. Furthermore, a 
quick service mode will be available at higher traveling cost to the 
courier service personnel. 

Application 5.6 : Sequencing of a drilling job 

A workpiece requires several holes of different diameters to be 
drilled on it. After drilling of each hole, we have to reset the machine 
(i.e. changing and positioning of the drill tool, clean the workpiece 
etc.), incurring a set up time. Let C denote the set up time between 
drilling of two holes i and j. The objective is to find out the optimal 
processing sequence of the holes that will minimize the total set up time. 
The problem can be formulated as the traveling salesman problem - the 
"machine" which functions as the "salesman", needs to visit the jobs in the 
time effective manner. This application will be similar to Application 5.4 
except constraint (5.4e) will not be present and the objective will be to 
minimize the total set up time. 

Application 5.7 : Interconnection of local area networks 
Reference : Fetterolf and Anandlingam (1991) 

Local area networks (LANs) are popular means of interconnecting 
information systems in corporate environments. LANs are used to 
interconnect PCs, workstations, mini and mainframe computers, fileservers, 
database servers and printers to share the benefits of multi-user systems. 

In big organizations, multiple LANs (as many as 10 or 20) of varying types 
(Ethernets or token rings), capacities and costs are connected through high 
speed devices, called bridges, to monitor the transmission of datapackets 
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among the LANs. 

The problem of optimal interconnection of LANs deals with the design 
of a LAN-LAN internetwork, with a given set of LANs, a traffic matrix 
defining host to host internetwork traffic, and a set of bridges of varying 
cost and capacities that will minimize the total cost of wiring and bridges 
costs. It will also satisfy following requirements : 

a) Network flows should satisfy traffic matrix, i.e., requirement of each 
user to receive or send data must be satisfied. 

b) Total average flow of data packets into each LAN should be below the 
recommended maximum for that LAN and it should also be below the 
maximum filter rate of each bridge connected to that LAN, i.e., 
maximum capacity of each bridge to determine whether packets should be 
transferred or not. 

c) The total average flow of data packets into each bridge should be 
below the maximum transfer capacity of that bridge between two LANs. 

d) To prevent data packets from getting caught in loops, internetwork of 
LANs should be a spanning tree. 


Let L denote the set of LANs and K denote the set of bridges of 
different capacity and cost. LANs are modeled as nodes of networks, and 
bridges are modeled as a set of arcs [ (i, j) , ( j, i) ] . As the LANs are fixed, 
node capacities are fixed. However, arc capacity is dependent on the type 


of bridge used between two LANs. Let t_ represents inter LAN traffic from 
LAN i to LAN j. Furthermore, a datapacket that is sent from LAN o to LAN d, 


is defined as a commodity od and represents the average flow of 

commodity od over the bridge connecting LAN i and LAN j, such that the 


total flow across the bridge (i,j) will be : 



i i 

oeL d€L 


x od 

ij 



58 


Let (Lj be the wiring and installation cost associated with using a 

k 

bridge to connect LANs i and j and C be the cost of purchasing a bridge of 
capacity k. We define a zero-one variable Y_ to denote whether a bridge 
of capacity k is placed between LANs i and j or not. 


The problem will be as follows : 


min E Z E 

keK ieL j>i 


k k 

(C. . + C K ) Y . . 
1J ij 


subject to 


Z 

jeL 



Z 

jeL 



t , if i = d (destination) 
od 

■ -t if i = 0 (origin) 

^ 0 otherwise 


for all i e L, for all o e L, for all d e L, 

(5.7a) 

X + X ^ I g k Y k for all i e L, j > i, (5.7b) 

J J keK J 

E (x. . + t .. ) i (£ k - h.) I(Y k ) + h. for all jeL, 
ieL 1J J1 J jn J 

for all n e L, for all keK, (5.7c) 

Z Z Z Y k = | L | - 1, (5. 7d) 

keK ieL j>i J 

X. . = T T X?^ for all ieL, for all j e L, (5.7e) 

ij it L jiL ij 

X?^ i 0 for all ieL, for all jeL, for all o e L, for all d e L, 
ij 

(5. 7f ) 

Y k j e {0,1} for all ieL, for all jeL, for all keK. (5.7g) 
Constraint (5.7a) are network flow conservation constraints. Constraint 

k 

(5.7b) restricts the total flow across each bridge, where g is the bridge 
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transfer rate. In constraint (5.7c), l is the filter rate of bridge k and 
hj specifies the capacity of LAN j. This constraint restricts total flow 
in each LAN. Constraint (5.7d) ensures the spanning tree structure of the 
inter connection network. 

If we relax (5.7b) and (5.7c) by using two Lagrangian multipliers A. 

2 

(for constraint (5.7b)) and A. , (for constraint (5.7c)), the relaxed 

jnk 

problem can be decomposed in two subproblems : 

Subproblem (1) : 

min E E E < c ij ♦ A V*. - I E e xj g k A 

keK ieL j>i J J ieL j>i keK J J 

- E E E *1* « k - h ) i(v k ) 

jeL neL keK J J J 

subject to (5.7d) and (5.7g). 

To use the minimum spanning tree algorithms, it is essential to 

determine the cost of connecting i and j relative to subproblem 1. Let A_ 
denote this cost for any arbitrary pair of LAN i and LAN j. The spanning 
tree requirement dictates that only one bridge type in the set K may be 
used to span i and j. Therefore the minimum cost A_ will be 


A. . = min -f C. . + C k - A 1 . g k - A 2 (l k - h. ) - A 2 (l k -h J 
1 J keK \ ij iJ ijk i Jik j 


This subproblem will be a spanning tree problem. 

Subproblem (2) : 

min E E X? (X . X ) + E E E E X L «ij * tji*- 

ieL j>i 1J 1J J1 jeL neL keK ieL JnK 1J J1 
subject to (5.7a), (5.7e) and (5.7f). 

This will be an uncapacitated multicommodity flow problem. This can 


also be viewed as minimum cost flow problem for each commodity od. In 
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order to minimize the cost of flows between 0-D pairs on an uncapacitated 
network, it is required to route each 0-D flow on the minimum cost path 
(shortest path) between origin and destination. Let denote the 
distance (cost) between each pair of nodes, such that 


*!. 

+ E 

E 

A 2 v 

if i < j. 


neL 

keK 

jnk 


A 1 . 

+ E 

E 

A 2 V 

if j < i. 

■ Ji 

neL 

keK 

jnk 



The shortest path between all pairs of nodes can be calculated by using all 
pair shortest path algorithms. 

Application 5.8 : Electrical distribution system design 
Reference : Current and Pirkul (1991) 

A national electricity board has to determine the design of an 
electrical distribution system. They want to connect two major power 
generating plants with high voltage transmission lines to enhance the 
reliability of the system and to facilitate load sharing. Furthermore, to 
connect the end users of electricity with this transmission system, lower 
voltage distribution lines are required. Transformers and electrical 
substations at the interconnection of high voltage transmission line and 
lower voltage distribution lines are also needed. The problem of optimal 
electrical distribution system design identifies a two-level hierarchical 
network that includes a path of primary arcs (higher voltage line) between 
two power generating plants and all the demand nodes for electricity must 
either be at the site of a transshipment facility (transformer etc.), or to 
be connected to such a facility via a path of secondary arcs (lower voltage 
lines). All transshipment facilities must be located on the primary path. 

We define an undirected network G = (N, A) with N as a set of n nodes, 
where node 1 and node n represent the power plants and other nodes 



61 


represent the demand nodes for the electricity, and A as a set of m arcs. 


Let CLj and C_ be the primary and secondary arc costs of each arc (i,j)€A. 
Furthermore, fj denotes the cost of constructing a transshipment facility 
at location j. 


We define a zero-one variable j indicating whether a primary arc 

connects node i to node j and also another zero-one variable Y. . 

J ij 

indicating whether a secondary arc connects node i to node j or not. We 
also add one dummy node n+1 in the network G(N,A) such that the new network 
will be G 7 = (N 7 ,A 7 ), where 


N 7 = N u n+1, 

A 7 = A u A, and A = {arc (i,n+l)|i e N}. 


We define arc cost C 7 „ = f. for 

i,n+l 1 

formulated as follows : 

n-1 n n n+1 

min Y Y C. . X. . + Y Y 

Till *- J 

j = 2 J J i=i j=i 
subject to the constraints : 


n 




E 

j=2 

x u 

= 1 , 


n-1 




E 

i=l 

X. 

in 

= 1 , 


n-1 


n 


E 

i=l 

X. . 
ij 

- E x jk = 0. 

k=2 JK 

for j 

E 

ieQ 

E 

jeQ 

X.j =5 |Q| - 1 

for all 

n+1 




E 

j=i 

Y. . 
ij 

= 1. i = 1,... 

,n. 

j*i 




n+1 


n 


E 

j=i 

Z. . 
ij 

- E z ki = i. 

k=i K1 

i = 1, 

j*i 


k*i 



all i e N. The problem can be 


(5.8a) 

(5.8b) 

= 2, . . . ,n-l, (5.8c) 

Q Q N such that j Q J ^ 2, (5.8d) 

(5. 8e) 

...n, (5 . 8f ) 
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Z.j - (n) Y s 0 i = 1 ,n, j = 1 n, j * i, (5.8g) 

n-1 

Y j,n+l " .E x ij " 0 J = 2 n-1. (5 . 8h) 

X^j e (0,1) for all (i,j) e A , (5.8i) 

Z ij . £ 0 for all (i,j) € A', (5.8j) 

Y.j e (0,1) for all (i,j) e A', (5.8k) 


{ 1, if a transshipment facility is located at i 
0, otherwise. 


Constraint sets (5.8a) - (5.8d) ensures that the solution will include a 
simple path of primary arcs from node 1 to node n. If dummy arcs = 
1, then a transshipment facility is located at node i; therefore, 
constraint (5.8h) requires the transshipment facility located on the 
primary path. Constraint sets (5.8e) - (5.8g), (5.8j) - (5.8k), ensures a 
spanning tree structure of secondary arcs on G' rooted at node 1. 

The problem can be solved by relaxing the constraint set (5.8h). We 
obtain Lagrangian subproblem as follows : 


n-1 n 


n n+1 


n 


n-1 


(L) min Y Y C. . X. . + Y Y C' . Y. . + Y X . (Y . - Y x. .). 

i=i jia l J ^ i=i j=i ^ 'J J J’ n+1 1=1 


subject to the constraints (5.8a) - (5.8g),and (5.8i - 5.8k). The 


Lagrangian subproblem can be decomposed into the following two subproblems: 


n-1 n 

(L ) min I l (C - A ) X . 

i=l j=2 J J J 

subject to the (5.8a) - (5.8d), (5.8i). 
n n+1 

(L ) min £ E C Y . . 

i=l j=l J J 
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subject to the (5.8e - 5.8g), (5.8j) - (5.8k), where 




f C' if (i,j) € A , 

I C'. . + X. if (i, j) e A' . 
1 ij J J 


Here ^ 0 is a Lagrangian multiplier. 

The problem has the structure of a shortest path problem and the 


problem has the structure of a minimum spanning tree problem. 


Application 5.9 : Design of mass transit system 
Reference : Current and Pirkul (1991) 


A state transport corporation wants to design the public mass transit 
system between two major cities. The primary mode might be a fixed rail 
system such as a subway, while the secondary mode might be a feeder system 
consisting of buses. Passenger stations as transshipment facilities are 
required at the locations on the primary subway where bus system will feed 
into subways. The problem of optimal design of mass transit system is also 
a two-level hierarchical network design problem like Application 5.8. We 
have to design the network such that there should be one primary rail 
system between two major cities and other cities should be connected either 
via a feeder bus system or should lie on the primary path. 

Application 5.10 : Optimal design of oil/gas pipeline network 
Reference : Current and Pirkul (1991) 

Oil and gas commission of a country wants to connect a major oil/gas 
field with a transporting facility or a processing facility. Furthermore, 
in that geographical area several other smaller oil/gas fields are also 
present. The OGC wants to design a two-level hierarchical network of 
pipelines at minimum cost to transport all available oil/gas to the 
processing centre. Typically, a main pipe line connects the major oil/gas 
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field with the processing centre. Small oil/gas fields in the area are 
connected to the main line via secondary lines, which are called gathering 
lines and consist of pipes having much smaller diameters. The problem can 
also be formulated like Application 5.8. 

Application 5.11 : Optimal transportation network design for resource 

extraction area 

A big industrial house is planning to set-up a new processing centre 
to process Boxite for manufacturing aluminum. They are planning to set-up 
this plant near the boxite extraction area which consists of scattered 
boxite mines. The company wants to connect the processing centre to the 
national transport system, with a nearest site on the national highway to 
facilitate transportation of finished products. Furthermore, all mines in 
that geographical area have to be connected with processing centre to 
facilitate transport of raw material. 

The problem is also a hierarchical network design problem with a 
primary path of railroad between processing centre and nearest site on 
national highway and each mine must either be on this primary path or be 
connected to primary path via a secondary transportation system (say 
trucks). Loading and unloading facilities are required at the 

interconnection of primary and secondary arcs. The problem will be 
similar to that in Application 5.8. 

Application 5.12 : Optimal design of transportation network of a 

developing country 

Reference : Current, Re-velle and Cohon (1986) 

In a developing country, transportation network between two major 
cities has to be designed. The primary path could be thought as an all- 
weather highway. The secondary paths might be unimproved roads connecting 
intermediate villages to the all-weather highway. All intermediate village 
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allocated among various activities and (for k = 1,...,K) denotes the 
total availability of resource k. Let d^ denotes the duration of activity 

i. The problem can be formulated as follows : 
min Y t Y . 

t nt 

subject to 


E Y i t 1 ’ 1 1 ’ ’ ' ' ’ n * 


I t(Y jt ~ Y it } d i € A > 


(5.14a) 

(5.14b) 


I r, 


ik 


( ^ 

H-d.+l 


Y. 

im 


^ b 


k ’ 


k = 1 , . . . ,K, 

t = 1, . . . ,T, 


Y it e 


(5. 14c) 
(5. 14d) 


Constraint (5.14a) ensures that every activity must start once. Constraint 
(5.14b) are the precedence constraint, where A is the set of all activities 
with precedence constraints and constraint (5.14c) are the resource 
constraints. Relaxing the resource constraints (5.14c) of the above 
formulation with Lagrangian multipliers X £ 0, the relaxed problem will 
be : 


min £ 
t 


t Y 


nt 


E Z "it v u 

1 t 


- E E 

k t 


A tk V 


subject to 


where 


EV U = 1, i = 1 n, 

E t( Y jt ‘ V ed i ( 1 'J> * A ' 

Y it * 


t.-d +1 
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The relaxed problem can be viewed as a generalization of a longest path 
computation in which we have to minimize not only the completion time but 
the sum of ‘costs’ of starting activity i at time t. Christofides 

(1987) suggested an algorithm to solve the relaxed problem. 

Application 5.15 ; Procurement of aviation fuels 
Reference : Austin and Hogan (1976) 

Defense department requires aviations fuels for several military 
installations. Generally, requirements for aviation fuels are met with 
purchases made in the usual competitive bidding environment. Invitation 
for bids are submitted to oil companies at frequent intervals annually. 
The objective of defense department is to procure the required fuel at the 
lowest laid-down (fuel plus transportation to demand sites) to the 
government. Generally, a wide variety of bidding options are available to 
the oil companies, such as availability of different transportation modes 
(pipe-line, tank car) or an oil company may indicate a total maximum offer 
which is less than the sum of maximum available at Its shipping points. 
Furthermore, a company may "tie in" its bid to one or more separate 
companies, that is, it may make its maximum offer contingent upon the 
amount awarded to other bidders. Sometimes, companies may indicate a 
minimum acceptable quantity for one or more shipping points. 

The problem has been formulated as a resource constrained minimum cost 
circulation problem by Austin and Hogan [1976]. Because of the features of 
"tie-in" between two companies, the total flow on some of the arcs must not 
exceed a limited quantity. This feature was dealt as a side constraint in 
the network. The additional side constraint was relaxed using Lagrangian 
relaxation technique and the resulting subproblem was a minimum cost 


circulation problem. 
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Application 5. 16 : Optimal planning for a forest industry 
Ref erence : Jornsten and Vabrand (1976) 

A forest industry is having various plants at different locations 
where at sawmills, timber logs are converted into timber and other by 
products such as pulp and paper in given proportions. There are several 
customers for these products and timber, which are having different but 
known demands for each of these. Furthermore, timber logs are supplied 
from different forests having known availability of timber logs. Timber 
can be either transported directly to the customers from the forests or to 
a plant where the timber is converted into pulp and paper. The cost of 
transportation of timber from different forests to customers or to the 
industries is known. Furthermore, transportation cost of pulp and paper 
from any plant to customers is also known. 

The objective of optimal planning for a forest industry problem is to 
find out a minimum cost transportation plan of the timber and by-products, 
keeping in view the demands at various customer centres, supply of timber 
from different forests and conversion proportion of by-products at 
different plants. 

The problem was formulated as a multicommodity generalized flow 
problem by [Jornsten (1986)]. A Lagrangian relaxation of complicating 
conversion constraints at each plant resulted into a convex cost flow 
subproblem. For a complete study of the model, refer to Jornsten and 
Varbrand (1986). 

Application 5.17 : A logistics planning system 
Reference : Klingman, Mote and Philips (1986) 

A logistics planning system was developed for W.R. Grace company, one 
of the largest suppliers of phosphate-based chemical products of the USA by 
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[Klingman et al. (1986)]. The mathematical model underlying this system 
includes production, distribution policies of multicommodities in multiple 
time periods. They decompose the problem into a generalized network 
component and a small linear non-network component. The generalized 
network component was transformed into pure network and non-network 
constraints were relaxed by using Lagrangian relaxation technique. The 
relaxed problem was solved as a pure network flow problem. As the model is 
tidious and complex, it is difficult to present it in our report, we 
decided to omit the description of the model. For a complete study of the 
model, refer to Klingman et al. [1986]. 
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